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INTRODUCTION 

The  techniques  for  analytically 
describing  the  pressure  gradient 
which  exists  from  the  static 
conditions  in  the  breech  to  the 
turbulent  area  directly  behind  the 
projectile  have  been  one  of  the 
central  problems  of  interior 
ballistics  from  the  lnceptiOTj  of  the 
science,  to  the  present  day.  The 
closed  form  interior  ballistics  of 
the  1940's  generally  used  a  constant 
factor  pressure  decrement,  based  on 
the  propellant  charge  to  projectile 
mass  ratio.  These  factors  were 
derived  from  an  integral  approach 
and  applied  incrementally.  Although 
they  were  nowhere  correct  on  a 
stepwise  basis,  in  the-  aggregate 
they  were  effective  and  suitable  for 
the  closed  form  calculations  of  the 
time. 

In  England  the  approach  used  by 
Corner  [Ref.  1]  and  Hunt-Hinds  [Ref. 
2]  was  that  of  a  shot  base  pressure 
P8  to  breech  pressure  ratio  Pb 
relationship  of 


ps  -  Pb/M  ♦  V2c/m) 

With  a  C/M  -  .5  this  predicts 
ps  “  -8pb 

In  the  United  States,  the 
Hirschf elder  approach  [Ref.  3]  used 
an  average  pressure  Pave  expression 
of  the  form 

P  .  P  (  3  ) 

rave  s  v  m 

which  with  a  C/M  -  0.5  predicts 
P_  -  .87  Pave  or  in  terms  of  Pb 
ps  -  -76  Pb 

Both  approaches  recognized  that 
these  approximations  became 
increasingly  less  accur  te  as  the 
charge  to  mass  ratio  increased  above 
0.5. 

As  computers  became  generally 
available  in  the  early  1960's,  other 
approaches  and  relationships  were 
p  stulated  for  the  lumped  parameter 
interior  ballistic  formulations 
which  were  developed  to  mod°l  the 
interior  ballistic  process  while 
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avoiding  solving  coupled  partial 
differential  equations.  An  approach 
used  for  many  years  by  one  of  the 
authors  [Ref.  M]  postulated  that  the 
pressure  differential  from  the 
breech  to  the  base  of  the  projectile 
could  be  expressed  in  terms  of  the 
velocity  of  the  projectile,  by  means 
of  "quasi-isentropic"  approaches, 
that  is: 


2  -y/y-i 

JL  ri  tH,  rLi 
-  [1  +  (  2  )  <YgRT)] 


or 

P 

s 


U 


-Y/Y-1 


.  [1  «■  (— )  (-K — )] 
Pave  r«RT 


with  the  value  Up/YgRT  being  the 

Mach  number  squared  of  the 
propellant  combustion  gases  at  the 
base  of  the  projectile. 

With  a  C/M  -  .5  this  approach, 
when  integrated,  would  give  an 

overall  effective  P  /P.  ratio  in  the 
.75  to  .8  range  similar  to  the  old 
Corner  and  Hirschfelder  approach. 

However,  as  C/M  values  of  very 
high  performance  gun  systems  started 
to  approach  and  even  exceed  values 
of  C/M  -  1  the  velocity  pressure 
postulate  was  found  to  theoretically 
under  predict  experimental  gun 
firing  experience. 

More  specifically,  as  in-bore 
velocity  Doppler  radar  systems 
became  available  in  the  early  1970's 
and  laser  based  visar  in-bore 
velocity  systems  in  the  mid  1970's, 
it  was  experimentally  seen  that 
substantially  more  AV  was  imparted 
to  the  projectile  later  in  the 
ballistic  cycle  than  predicted  by 
the  velocity  based  "quasi- 
isentropic"  postulate. 

Clearly  a  more  accurate 
postulate  was  required.  From  1 985 
to  1987  the  Air  Force  Armament 


Laboratory  sponsored  and  the  Air 
Force  Office  of  Scientific  Research 
funded  an  effort  at  the  University 
of  Florida  to  study  high  velocity 
in-bore  gas  dynamics  and  generate  a 
better  postulate  as  to  the 
mechanistic  description  of  in-bore 
gas  gradients. 

The  balance  of  this  paper 
provides  the  results  of  the  study 
and  indicates  that  the  in-bore 
pressure  gradient  is  not  driven  by 
the  velocity  of  the  projectile  but 
rather  by  the  acceleration  of  the 
projectile. 

PHYSICAL  MECHANISMS  CAUSING 
PRESSURE  AND  DENSITY  GRADIENTS 

The  physical  mechanism 

postulated  to  cause  a  pressure 
difference  between  the  breech  and 
projectile  is  the  acceleration  of 
the  projectile.  As  the  projectile 
accelerates,  it  produces  a  series  of 
expansion  waves  behind  it.  The 
pressure  drops  from  breech  to 
projectile  through  these  waves. 

The  ballistic  cycle  begins  as 
combustion  occurs  within  the  chamber 
bounded  by  the  breech,  tube  wall, 
and  projectile  and  increases  the 
pressure.  Because  the  burning  rate 
of  the  propellant  increases  with 
increased  pressure,  high  pressures 
are  quickly  produced.  These  high 
pressures  accelerate  the  projectile 
producing  a  series  of  exansion  waves 
behind  it.  These  wave  accelerate 
the  gases  behind  the  projectile  to 
the  velocity  of  the  projectile. 
They  travel  toward  the  breech, 
reflect,  and  return  to  the 
projectile  distributing  the 

velocities  of  the  gases  between  the 
breech  and  the  projectile. 

As  the  projectile  accelerates, 
the  gases  behind  the  projectile  also 
accelerate.  This  acceleration 

produces  a  gradient  of  pressure, 
density,  and  temperature  behind  the 
projectile.  Subsequently,  the 
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pressure  at  the  base  of  the  projec¬ 
tile  is  lower  than  that  at  the 
breech. 

The  pressure  drop  is  then 
governed  by  the  acceleration  of  the 
piston.  After  the  acceleration 
stops,  the  pressure  drop  from  the 
stationary  end  of  the  cylinder  to 
the  moving  piston  goes  to  zero.  In 
a  typical  ballistic  system  the 
acceleration  begins  low,  it  then 
increases  rapidly  to  a  peak  and 
drops  off  with  movement  of  the 
projectile.  It  would  then  be 
reasonable  that  the  drop  of  pressure 
from  breech  to  projectile  would  be 
dependent  on  the  acceleration 
history  of  the  recent  past.  The 
recent  past  would  be  defined  as  the 
time  between  formation  of  an 
expansion  wave  and  when  that  wave 
has  reflected  and  returned  to  the 
piston.  Once  the  expansion  wave  has 
returned  to  the  piston,  its 
remaining  effect  is  only  to  lower 
the  pressure  of  the  gas  in  the 
chamber . 

EQUATIONS  OF  MOTION 

FOR  UNSTEADY  FLOWS 


The  breech  to  projectile 
gradients  of  pressure,  density  and 
temperature  are  gorverned  by  the 
conservation  equations  of  mass, 
momentum  and  energy.  One 
dimensional  forms  of  these  equations 
were  solved  using  simplifying 
assumptions  to  determine  analytical 
relationships  between  the  projectile 
acceleration  and  the  gradients.  The 
equations  of  motion  for  a  mixture  of 
gas  and  solids  are  the  conservation 
of  mass  and  the  conservation  of 
momentum. 


at  (psi  °“v,)  +  py,)  + 


3x  ^psl^1-v'^  US1  +  Pv’u)  -  o  (1) 
ft  (usl  psl  (W)  +  upy,)  + 
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u  pv' 


It  was  assumed  that  the  velocity 
of  the  solids  is  related  to  the 
velocity  of  the  gas  by  a  constant 
ratio  throughout  the  tube  at  any 
particular  instant.  No  friction  or 
heat  transfer  effects  were 
included.  The  two  equations  were 
combined  to  give  an  equation  which 
has  a  form  similar  to  the  unsteady 
Bernoulli  equation. 
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The  terms  on  the  left  hand  side 
of  Eqn.  3  are  acceleration  terms. 
The  integrated  acceleration  produces 
the  pressure  difference  from  breech 
to  projectile.  The  right  hand  side 
represents  the  pressure  gradient 
produced  by  the  acceleration. 

Acceleration  Term 


The  acceleration  term  was 
evaluated  for  solids  and  gases.  A 
linear  profile  of  velocity  with 
distance  was  assumed  for  both  the 
velocities  of  gas  and  solids.  The 
velocity  of  the  solids  was  assumed 
to  lag  that  of  the  gas  by  the 
ratio  u^^/u.  This  approach  was 
used  because  several  different 
models  including  the  XNOVA  code 
[Ref.  5]  indicated  that  the  gas 
velocity  profile  was  quite  linear. 

When  the  acceleration  of  the 
piston  is  zero,  the  linear  velocity 
model  seems  to  predict  that  the 
velocity  past  a  point  is  always 
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decreasing.  Vlhen  a  >v2/L,  the 
velocity  past  a  point  will  increase 
with  time.  When  ap<v2/L,  the 
velocity  past  a  point  will  decrease 
with  time.  This  does  not  mean 
however,  that  the  velocity  of  a 
particle  is  slowing  down.  As  one 
follows  a  particle,  it  continues  on 
at  the  same  or  greater  velocity 
depending  on  the  level  of 
acceleration  of  the  piston. 

Assuming  a  linear  velocity 
profile  and  integrating  the 
acceleration  side  of  equation  3  gave 


p  ,  L  a  u 
r  ml  P  .  _P 

Lp  2  2 


p 

f  t  dP 
Jp  n 


afai,2  .  (fsi)m2 

pm  pm 


Pb  pm 


The  most  significant  portion  of  the 
acceleration  side  is  L  a ,12.  The 
factor  Pml/Pm  accounts  for  solids 
which  are  not  accelerating  as  fast 
as  the  gas.  The  projectile  velocity 
squared  portion  of  the  acceleration 
side  will  always  have  a  minus  sign 
in  front  of  it  because  p  ,/p  is 
always  less  than  1 .  When  there  are 
no  longer  any  solids,  p  ./p  goes  to 
1  and  the  entire  acceferafton  term 
goes  to  L  a_/2. 

The  ratio  p  1/p  is  determined 
knowing  the  pohos'lty,  the  gas 
density,  and  ugl/u. 

Ugi/u  was  estimated  as  a  constant 
along  the  barrel  at  any  instant  in 
time.  It  was  assumed  to  change  with 
time.  The  ratio  of  average 
acceleration  of  the  solids  over  the 
average  acceleration  of  the  gas  was 
assumed  to  be  equal  to  ugl/u.  The 
drag  force  on  the  cylindrical  grain 
was  assumed  proportional  to  (u-ugl)2 
and  equal  to  the  acceleration  of' the 
grain  times  its  mass. 


usl/u  was  solved  fo  r  as  a 
function  of  a  nondimensional 
acceleration  number.  Fa.  Fa  has 
been  formulated  to  be  a  function  of 
the  easily  determined  parameters. 


2CppU 

3lirop3i  aP 


1  +  2Fa  (1  ■  /4  Fa  +  1)  (6) 


The  assumption  of  a  linear 
velocity  profile  was  tested  using  an 
inviscid  computer  model  of  gas  flow 
in  a  tube  between  breech  and 
projectile.  The  projectile  was 
driven  by  estimates  of  accelerations 
from  experimental  shots  536  and  571 
which  will  be  discussed  later.  The 
ratio  of  the  integrated  particle 
acceleration  (left  hand  side  of 
equation  3)  divided  by  a  L/2  was 
plotted  against  time  in  Tigures  1 
and  2.  If  the  velocity  profile  was 
linear,  the  ratio  would  be  1  .  In 
Figure  1  the  ratio  is  0.98  during 
the  last  70$  of  the  time.  This 
indicates  that  the  integrated 

acceleration  estimate  of  a  L/2  is 
2$  high.  Figure  2  is  for  a  case 
with  a  30$  higher  peak  pressure. 
This  higher  pressure  causes 
increased  acceleration.  This 

increased  acceleration  rate  is  the 
probable  cause  of  the  hump  in  Figure 
2  that  starts  at  0.002  seconds. 
The  peak  of  the  hump,  1  .07, 

indicates  that  the  integrated 

acceleration  estimate  of  a  L/2  is 
7$  low.  P 

Pressure  Gradient  Term 

The  pressure  gradient  term  was 
evalated  for  the  two  cases  of  gas 
only  and  of  gas  and  solids.  The 
pressure  was  integrated  by  assuming 
an  isentropic  relationship  between 
pressure  and  density. 

The  integration  of  the  gas  only 
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case  stated  in  terms  of  breech 
temperature  and  density  directly 
yielded  equation  7a.  Y _ 

2  Y_1 

Pr  (Y-1 ) (1-p  n)  L  a  c 

—2  «  (i - 2 - E - ) 

PK  2Y  R  T  ' 

b  b 


When  stated  in  terms  of  breech 
pressure  and  density  it  yielded 


(r-1)  P„  La  E 


The  gas  and  solids  case  was 
integrated  using  a  linearized  form 

of  the  term  (P/Pb)^1/Y^.  The 

equation  was  statea  in  terms  of 
breech  temperature  and  density 

?F  p  (1/v’-1) 

•  1  ♦  Y  (-£±—- -  ♦  1)* 

b  pb 


[exp  ( 


-v*(1-pbn)Lap  C 


■)  -  1]  (8a) 


It  was  also  stated  in  terms  of 
breech  pressure  and  density. 

P  p  (1/v'-1) 

“  *  1  +  Y  (- -  +  1  )* 

b  pb 


[exp  ( 


-V  pbLap  t 

2YP. 
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)  -  i] 


The  multiplier  term 
(p  .  (1/v’  -  1)/pb  ♦  1)  can  greatly 
increase  the  pressure  drop  from 
breech  to  projectile  over  the  gas 
only  case.  This  is  because  the 
large  forces  are  needed  to 
accelerate  the  solids  which  have  a 


density  much  greater  than  the 
gases.  The  average  pressure  can  be 
determined  by  integration  of  a 
series  expansion  of  8b. 
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Kinetic  Energy  in  the  Gases 

The  ratio  of  kinetic  energy  in 
the  gas  solid  mixture  to  the  kinetic 
energy  of  the  same  mass  moving  at 
the  projectile  velocity  is  1/6.  6 
was  obtained  using  a  linear  velocity 
distribution  and  a  parabolic  density 
distribution  obtained  from  the 
pressure  distribution  of  Eqn.  8. 


p  (1-v')  p 

-  +  (1  - 

Pkv  3ph 

6  =  3  — - ~ - - -  (10) 

p  , (1-v' )  Uc  3p 
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The  case  of  gas  only 

6  3(1  ~  k 
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COMPARISONS  OF  ANALYTICAL 
MODEL  WITH  EXPERIMENTAL  DATA 


The  equations  just  described  in 
were  compared  with  experimental  data 


I  |.»  !.•  #«  |,»  V  tl‘  Ik*. 


from  Eglin  AFB  for  firings  with 
muzzle  velocities  of  1000  m/s. 
Comparisons  were  also  made  with 
computer  modelled  predictions 
obtained  from  Eglin  AFB  and  the 

Ballistics  Research  Lab.  The 
experimental  data  came  from  two 
series  of  shots  of  a  30  mm  aircraft 
cannon.  Each  series  consisted  of 

four  shots.  These  runs  were 
evaluated  by  shifting  the  time 
scales  to  make  the  peak  pressures 
coincide  and  then  comparing  the 

measured  pressures  in  time.  The  two 
pressure  time  traces  which  best 

matched  each  other  were  then  used  as 
the  standard  data  for  each  series. 

Each  set  of  data  consisted  of 
time  and  for  pressure  readings.  The 
first  reading  was  called  the  chamber 
pressure  reading.  It  was  not  taken 
at  the  breech  but  just  in  front  of 
the  initial  position  of  the  rotating 
band.  The  three  other  readings  were 
taken  at  three  downbore  locations. 

The  data  in  a  nondimensional 
form  were  used  for  the  comparison 
with  analytical  equations  (7  and 
8).  The  pressures  at  the  three 
downbore  locations  were  normalized 
by  the  chamber  pressure  reading. 
The  resulting  pressure  ratio  was 
less  than  1  because  the  downstream 
pressures  were  always  lower  than  the 
chamber  pressure. 

The  analytical  pressure  ratio 
was  determined  using  equation  8  to 
find  the  downbore  pressure  with 
respect  to  the  breech  pressure  for 
both  the  chamber  location  and  the 
downbore  location.  Values  of  the 
variables  in  equation  8  were 
obtained  from  an  interior  ballistics 
program  of  Heiney  [Ref  4].  The 
variables  needed  were  average 
porosity,  breech  density  and 
pressure,  projectile  acceleration, 
and  breech  to  projectile  length. 
The  analytical  pressure  ratio  was 
determined  as  a  function  of  time  for 
each  downbore  pressure  tap  location 
beginning  from  the  time  that  each 


tap  was  passed  by  the  projectile. 

The  results  are  quite 
encouraging.  Comparisons  with  the 
536  experimental  data  is  shown  in 
Figures  3~5.  The  experimental  data 
is  the  jagged  line.  The  other  two 
lines  are  the  results  of  calcula¬ 
tions  from  equation  8  with  the  upper 
line  being  a  gas  only  calculation 
and  the  lower  line  being  the 
calculation  when  solids  are  includ¬ 
ed.  At  1700  microseconds  in  Figure 
3  one  can  see  that  the  solids  theory 
agrees  quite  well  with  the  data. 

At  a  later  time  of  2500 
microseconds  in  Figure  4  the  gas 
only  theory  seems  to  be  a  better 
estimate  than  that  of  the  solids. 
This  is  probably  due  to  the  fact 
that  the  solids  do  not  follow 
directly  behind  the  projectile  but 
lag  it.  The  £  squared  term  in 
equation  8,  causes  the  majority  of 
the  pressure  drop  to  occur  in  the 
30{  of  the  length  just  behind  the 
projectile.  If  the  solids  concen¬ 
tration  is  very  low  in  this  region, 
the  pressure  will  drop  less  and  will 
better  follow  the  gas  only  theory. 

The  last  figure  in  the  536 
series  is  Figure  5.  It  shows  good 
agreement  between  theory  and 
experiment  for  either  gas  only  or 
gas  and  solids.  The  two  theoretical 
lines  are  very  close  because  at  3000 
microseconds  there  are  very  few 
solids  left. 

Comparisons  of  571  experimental 
data  with  theory  are  shown  in 
Figures  6-8.  This  shot  had  a  muzzle 
velocity  7%  higher  than  the  previous 
536  shot.  This  was  produced  by  a 
thinner  grain  web  thickness.  Figure 
6  shows  the  pressure  ratio  for  down¬ 
bore  tap  2.  Burning  has  occurred 
more  quickly  and  few  solids  remain 
so  that  the  gas  and  solids  lines  are 
very  close.  Predicted  pressures  are 
lower  than  those  measured.  One 
explanation  may  be  lack  of  time  for 
the  velocity  profile  to  become 
linear  after  the  time  of  peak 


pressure. 

Figure  7  shows  quite  good 
agreement  with  the  experimental  data 
over  the  last  half  of  the  time  of 
the  run.  Figure  8  shows  a  *  'essure 
ratio  higher  than  the  data.  This 
may  be  related  to  the  hump  in 
acceleration  ratio  shown  in  Figure  2 

Two  higher  speed  cases  were 
compared  with  theory.  Experimental 
data  was  not  available  for  either  of 

these.  However,  computer  predic¬ 
tions  were  available  using  two 
different  codes  which  model  the  flow 
of  both  solid  and  gas  as  the  solids 
burn. 

The  first  case  had  a  1400  m/s 
muzzle  velocity.  Combustion  was 
complete  and  gas  only  theory  was 
compared.  Pressure  at  the 

projectile  base  over  pressure  at  the 
breech  was  compared  when  the 
projectile  wa3  at  the  muzzle.  The 
computer  prediction  of  the  pressure 
ratio  was  0,707  and  the  pressure 
ratio  of  the  theory  of  this  work  was 
0.665. 

The  second  case  was  a  large 
diameter  round  with  a  1450  m/s 
muzzle  velocity.  Combustion  was 
complete  and  the  pressure  ratios  was 
compared  when  the  projectile  was  at 
the  muzzle.  The  computer  prediction 
of  the  pressure  ratio  was  0.731  and 
the  pressure  ratio  of  the  theory  of 
this  work  was  0.747. 

In  both  of  these  cases  the 
theory  matched  the  computer 
predictions  well.  The  Mach  number  of 
the  gases  behind  the  projectile  in 
both  of  the  last  two  cases  is  near 
1.5  indicating  that  the  flow  was 
supersonic. 

SUMMARY 

The  gradient  in  pressure  and 
density  from  breech  to  projectile 
was  postulated  to  be  caused  by  the 
acceleration  of  the  gases  behind  the 
projectile.  The  major  governing 


parameter  was  the  projectile 
acceleration  times  the  axial 
position  of  the  projectile. 
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NOMENCLATURE 


ap 

acceleration  of 
projectile 

V 

- 

adjusted  acceleration 

cp 

- 

change  mass 

Fa 

- 

acceleration  number 

L 

—■ 

distance  x  between 
breech  and  projectile 

M 

- 

projectile  mass 

P 

- 

pressure 

p 

ave 

" 

average  chamber 
pressure 

R 

- 

gas  constant 

t 

- 

time 

T 

- 

temperature 

u 

- 

velocity  of  gas 

up 

- 

velocity  of  projectile 

kinetic  energy  of  mass 
of  gas  at  Up/kinetic 
energy  of  gasF 

y 

' 

specific  heat  ratio 
Cp/Cv 

n 

- 

covolume 

v’ 

- 

porosity 

6 

- 

x/L 

pm 

Pm1 

- 

average  density 

density  times  velocity 
average 

pm2 

density  times  velocity 
squared  average 

propellant  grain 
density 

psl’ 

- 

density  of  a  grain 

including  gas  within 
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Figure  2.  Instantaneous  Particle  Acceleration  Integrated 
from  Breech  to  Projectile  Over  L  ap/2  Calculated  Using 
Inviscid  Theory  of  a  Compressible  Gas  Using  Acceleration 
of  the  Projectile  from  Shot  571. 
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Figure  3.  Pressure  at  Dournbore  Tap  #2  Over  Pressure  at 
Chamber  Tap  fl  as  a  Function  of  Time. 


1.2 


1.00 


0.0' 
0.000 


0.001  0.002  0.003 

TIME  (SECONDS) 


0.004 


Figure  1.  Instantaneous  Particle  Acceleration  Integrated 
from  Breech  to  Projectile  Over  L  ap/2  Calculated  Using 
Inviscid  Theory  of  a  Compressible  Gas  Using  Acceleration 
of  the  Projectile  from  Shot  536. 
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Figure  4.  Pressure  at  Dovrnbore  Tap  #3  Over  Pressure  al 
Chamber  Tap  #1  as  a  Function  ol  Time 
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Figure  5.  Pressure  at  Downbore  Tap  #4  Over  Pressure  at 
Chamber  Tap  #1  as  a  Function  of  Time. 
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Figure  8.  Pressure  at  Downbore  Tap  #4  Over  Piessure  at 
Chamber  Tap  #1  as  a  Function  ot  Time. 
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Figure  6.  Pressure  at  Downbore  Tap  #2  Over  Pressure  at 
Chamber  Tap  #1  as  a  Function  of  Time. 
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Figure  7.  Pressure  at  Downbore  Tap  43  Over  Pressure  at 
Chamber  Tap  #1  as  a  Function  of  Time 


